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ABSTRACT 

We have used the SHARP polarimeter at the Caltech Submillimeter Observatory to map the polar- 
ization at wavelengths of 350 and 450 (im in a ~ 2' x 3' region of the Orion Molecular Cloud. The map 
covers the brightest region of the OMC-1 ridge including the Kleinmann-Low (KL) nebula and the 
submillimeter source Orion-south. The ratio of 450-to-350 /an polarization is ~ 1.3 ± 0.3 in the outer 
parts of the cloud and drops by a factor of 2 towards KL. The outer cloud ratio is consistent with 
measurements in other clouds at similar wavelengths and confirms previous measurements placing the 
minimum of the polarization ratio in dusty molecular clouds at A ~ 350 /im. 

Subject headings: dust, extinction — ISM: clouds — ISM: individual (OMC-1) — polarization - 
submillimeter 



1. INTRODUCTION 

Studies of the wavelength dependence of interstellar 
polarization are com mon place in the near- visible region 
of the spectrum (e.g.. | Serkow ski 1958; M artin fe WhitteD 
[19901: iMartin. Clavton. fc Wolfl 119991 : iWhittetJ I2004T ). 
Such studies have put con straints on the mechanisms of 
magn etic grain alignment fRobcrge 2004; Lazarianll2003L 
12007ft as well as the properties of dust grains responsible 
for interstellar extinction (e.g.. [Hildebrand fc Draeovanl 
ll995HAndersson fc Potterll2007tlWhittet et al.ll2008t and 
references therein). There are fewer studies of the 
wavelength dependence at far-infrared and submillime- 
ter wavelengths (~ 0.1 - 1 mm) where dust grains are 
responsible for most of the observed emission in Galactic 
clouds. 

The first studies of multi- wavelength far-infrared/sub- 
millimeter polarimetry found unexpected results. Rather 
than a featureless spectrum due to a single population 
of dust grains, the polarization spectrum is observed to 
fall from 60 to 10 to 350 iim before ri s ing again to 
850 and 13 00 /jm (|Hildebrand et all 119991 : IVaillancourt] 
l2002t 120071 ) . This spectral structure is attributed to the 
existence of multiple dust grain populations whose po- 
larizability or alignment efficiency is correlated with ei- 
ther the grain temperature, the spectral dependence of 
the emissivity (the "spectral index" ) , or a combination 
of both. However, the existing data are quite sparse in 
terms of both wavelength coverage and the types of ob- 
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jects observed (bright Galactic clouds). Therefore, rather 
than pursue observations with the ability to test specific 
physical models, our immediate goal is to improve the 
empirical description of the spectrum by measuring the 
polarization at additional wavelengths between 100 and 
850 /im which bracket the observed minimum at 350 /im. 

To further this goal we have begun a campaign to mea- 
sure the polarization of Galactic clouds at both 350 and 
450 /im using SHARP, the SHARC-II polarimeter at 
the Caltech Submillimeter Observatory on the summit 
of Mauna Kea. While this is a fairly short wavelength 
baseline, such measurements will further constrain the 
location of the polarization minimum. For example, one 
can simply ask whether the minimum is less than, greater 
than, or approximately equal to 350 /im. This will place 
additional constraints on the possible range of polariza- 
tions, temperatures, an d emissivities of the const ituent 
grain populations (e.g.. iHildebrand fc Kirbvil2004f ). 

We have carried out the first set of observations for 
this project towards the Orion Molecular Cloud (OMC- 
1), a bright and well-studied region of massive star for- 
mation (e g..lHoude et al.ll2004tlJohnstone fc Ballvi ri999: 
iLis et al.l 1 19981) . In §2] below we review the SHARP 
instrument and the polarimetric/photometric observa- 
tions. Maps are presented in followed by a discussion 
of the polarization spectrum in Sj4] 

2. OBSERVATIONS AND DATA REDUCTION 

SHARP (|Li et al.l [20061 . [20081 : iNovak et alJl200l is a 
fore-optics module that installs onto the SHARC-II cam- 
era. Incident radiation is split into two orthogonally po- 
larized beams which are then imaged onto opposite ends 
of the 12 x 32 pixel SHARC-II bolometer array. The re- 
sult is a dual-polarization 12 x 12 pixel polarimeter with 
a 55" x 55" field of view (FOV). The polarization is 
modulated by stepping a half-wave plate (HWP) at the 
relative angles 0, 22.5, 45, and 67.5 degrees. 

Within each FOV and HWP position, standard pho- 
tometric beam-switching was performed with a chop 
throw of 5' and a chop position angle in the range 80° 
- 150° east of north. This is repeated at positions 
separated by 50" in right ascension and declination to 
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build maps larger than the FOV. For this work we 
have generated maps of size 2' x 3' encompassing the 2 
brightest submillimeter cores in OMC-1 (the Kleinmann- 
Lo w nebula, hereafter KL; and the s u bmillim eter source 
of iKeene. Hildcbra nd. fc Whitcombl (|1982l ) sometimes 
called Orion-south, hereafter KHW). Observations at 
350 and 450 /im were carried out on 2007 February 16 
and 2006 December 5, respectively. 

The generation of Stokes parameters, polarization am- 
plitudes, and position angles within each single set of 4 
HWP a ngles follow s the same general procedures out - 
lined by lPlatt et all (|1991f ) and lHildebrand et all (|2000l ). 
Maps of the linear Stokes parameters I, Q, and U are 
generated from the dithered and stepped array positions 
by in terpolating the data onto a fi ner grid and coadding 
(e.g.. lHoude &: Vaillancourtj|2007j ). 

The polarization data pr esented in have been 
corrected for posi tive bias (|Simmons fc Stewart! 119851 ; 
iVaillancourtJ 120061 ) and for measured polarization effi- 
ciencies of 93% and 98% at 350 and 450 /on, respectively. 
The data have also been corrected for instrument polar- 
izations of » 3-0.5% at 350 and < 0.2% at 450 
(|Li et alJl2008l : Vaillancourt et al., in preparation). The 
polarization position angle calibration, measured by illu- 
minating SHARP with an unpolarized source through a 
calibration grid, is accurate to within 2°. 

We obtained photometric data at 350 /im using 
SHARC-II in camera mode on 2007 August 11-12. Four 
raster scans (without chopping) centered on Ori IRC2 
were performed to cover a ~ 10' x 10' field. Flux maps 
were produced using the sharcsolve 9 utility, which 
models the source flux, sky background fluctuations, in- 
strument gains, and drifts in the instrument electronics. 
Absolute fluxes are calibrated with respect to the stan- 
dard sources L 1551 and Mars (using peak flux estimates 
of 45.2 and 3.80 x 10 3 Jy per 9" FWHM beam, respec- 
tively) 10 and have an uncertainty of ~ 20%. 

3. RESULTS 

Figure [T] shows the polarimetric and photometric maps 
in the ~ 2' x 3' region covered by the polarization obser- 
vations. The size of the gaussian smoothing kernel used 
to interpolate the 350 /im polarization data has been cho- 
sen so that the resulting resolution matches that of the 
450 jum polarization data (13" FWHM). Figures [TJa-c 
plot the 350 and 450 /im polarization results over a con- 
tour map of total flux measured at 350 /im. At each 
wavelength we plot vectors at an interval of 9'.'5 (Figs. 

EEM). 

Figure [T]a shows the polarization results with the 
length of the vector proportional to the measured po- 
larization amplitude and the position angle parallel to 
the polarization vector. The most obvious feature of this 
map is the decrease in the polarization amplitude (at 
both wavelengths) towards the KL intensity peak, and 
to a lesser extent towards KHW. This so-called "po- 
larization hole" effect in which the polarization drops 
towards intensity peaks has been observed in these and 
other sources by experiments at a range of far-infrared 

9 http: //www. submm. caltech. edu/$\sim$sharc/analysis/overview. 
No significant differences were seen between maps reP 
duced usin g sharcsolve and those reduced u sing CRUSH 
(http: //www. submm. caltech. edu/$\sim$sharc/crush/ ) 
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P(450) / P(350) 

Fig. 2. — Histogram of the 450/350 /im polarization ratio. All 
data shown here have been limited to only those points where P > 
3rr p at both wavelengths. Also shown are histograms where data 
points satisfy the additional criteria that the position angle rotates 
by less than 10° between the two wavelengths (|A<£| < 10°) and 
that the points are at least 20" away from the two submillimeter 
flux peaks KL and KHW. 



1 http: //www. submm . caltech. edu/$\sim$sharc/analysis/calibrat 



and submillimeter wavelengths (e.g.. ISchleunind Il998t 
Coppin et ail 120001 : IMatthews et all 1200 H IDotson et all 
2000. IH3). 

In Figure [TJ> we have rotated the polarization vectors 
by 90° to show the inferred magnetic field direction. The 
plotted -B-vectors are drawn with a constant length (not 
proportional to polarization amplitude) in order to more 
clearly observe the position angle differences across the 
map and between the two wavelengths. The position 
angles exhibit a clockwise rotation with increasing wave- 
length, differing by more than 25° around KL and the 
region immediately east (less than 1% of the data exhibit 
larger changes). This angular rotation is much smaller 
to the west and south of KL. The larger variation of the 
position angle with wavelength towards KL than along 
the ridge between KL and KHW is also clearly seen when 
the data are compa red with observations a t 850 /im (not 
shown, but see e.g.. IMatthews et"aT]|2003h . 

The color scale of Figure [T]c shows the ratio of 
the polarization amplitude at the two wavelengths, 
P(450)/P(350). Most points in the map exhibit a two- 
point polarization spectrum which increases with wave- 
length, P(450)/P(350) > 1 (Fig. [2]). The most obvious 
exceptions are towards KL and the region to its north- 
west. KHW does not exhibit a similar drop in the polar- 
ization ratio. 

4. DISCUSSION 

4.1. Minimum in the Polarization Spectrum 

Examining SHARP data at only points outside the KL 
and KHW intensity peaks (beyond a 20" radius) we find a 
median in the distribution of F(450)/P(350) = 1.3 ±0.3, 
where the uncertainty represents the standard deviation 
of those points (Fig. This data point is consistent 
with a simple interpolation of data from other Galactic 
clouds at wavelengths of 350 and 850 /im (Fig. [3|). (Note 
that Fig. [3] also includes data at 850 /im fro m the Galac- 
tic cloud W51 (IChrvsostomou et alj| 2002l. not available 
the time of the IVaillancourti (|2002f l compilation; the 
comparison betwee n the 850 /im SCUB A data and the 
350 /im He rtz data (|Dotson et al.ll2008l ) is performed as 
ienMiaed in lVaillancourtJ ( 2002| ).) From the results of Fig- 
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Fig. 3. — Far-infrared and submillimeter polarization spectrum, 
normalized at 350 fim. The 450/350 fim OMC-1 comparison from 
this work is shown as a solid triangle. The 850/350 fim com- 
pa rison in W51 (open square s) is calculated from 850 fim data 
in IChrv sostomou ct al. (20Q3) a nd 350 fim da t a in I Dotson et alj 
020081) . All other data are from IVaillancourtl [120021) . The solid 
curve is a 2-component dust model (see text). 

ure [31 and assuming that the polarization spectrum in 
OMC-1 is similar to that of other clouds (i.e., it con- 
tinues to rise at wavelengths greater than 450 /im), we 
estimate that the minimum in the spectrum cannot oc- 
cur at wavelengths much less than 350 /im. 

4.2. Polarization Ratio in Different Environments 

Given the low polarization ratio towards KL compared 
with the rest of the cloud, it is natural to ask if this may 
be due to some unique physical condition as compared 
with the remainder of the cloud. At the very least, it 
is obvious that KL is the region with the highest flux 
density in the cloud. To examine this trend, we plot the 
450/350 polarization ratio as a function of 350 /mi flux 
F 350 . Figure Q] plots the polarization data as individ- 
ual points (dots) and combined into equal-sized logarith- 
mic flux bins (diamonds). At low fluxes (P350 < 50% 
of the peak flux) the polarization ratio clusters about 
P(450)/P(350) w 1.2 - 1.5 with a large scatter (this dis- 
tribution is best shown in the histogram of Figure [2|). At 
larger fluxes the binned data show a sharp drop in the po- 
larization ratio. The flux level at which this drop occurs 
is indicated by the thick contours in Figure [T]c. With few 
exceptions all data above this flux level lie in KL. The 
binned data in Figure [4] drop below P(450)/P(350) = 1 
for P(350) > 60%. At this flux level all points are within 
~ 20" of KL. 

Polarized emission from dust grains all with the 
same temperature results in a nearly constant polariza- 
tion spectrum in the w avelength range 50 - 2000 /jm 
( Hild ebrand et al.| [l999) . A spectrum such as that in Fig- 
ure[3 or points in Figure|4] with P(450)/P(350) ^ 1, re- 
quires the existence of at least two dust-emission compo- 
nents (jHildebrand et al.lll999HVaillancourtll2002f ). Such 
a polarization spectrum is modeled in Figure [3l This 
model has a minimum at 330 /im and uses optically thin 
dust with temperatures T = 20 and 50 K, spectral in- 
dices (3=2 and 1, respectively, 11 and peak fluxes in the 
ratio Pcoid/Phot = 0.5. The cold component is unpolar- 
ized while the warm component is polarized. While a 

11 An inverse T-0 r elation is observed to be typical for dust in 
molecular clouds fe.g. JDupac et al.ll200ll . I2003T) . 
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Fig. 4. — Polarization ratio (450 / 350 (im) vs. the total 350 fim 
flux (percentage of peak flux). Dots indicate individual measure- 
ments within the maps of Figure [TJj while diamonds represent an 
average of the plotted quantity over logarithmic bins of 350 fim 
flux. Error bars on the polarization ratio represent the standard 
deviation of the mean within each bin. 

two-component dust model for molecular clouds is phys- 
ically unrealistic, it is sufficient to illustrate how a multi- 
component model can explain the empirical results. 

Even with only 2 components, small adjustments to the 
parameters allow for better agreement between the model 
and the specific data points of Figure [3J A spectral en- 
ergy distribution (SED) dominated by cold unpolarized 
dust (Pcoid/Phot 3> 1) can reproduce the ratio observed 
in the OMC-1 cloud envelope, P(450)/P(350) » 1.3. 
Such an SED will appear nearly-isothermal if the mea- 
surement uncertainties are large and/or it is not suffi- 
ciently sampled in w avelength space. Th is is consistent 
with the analysis of IVaillancourtl (|2002t ) where no sig- 
nificant evidence is seen for multiple dust temperatures 
outside of KL or the M42 H II region, but where the 
polarization spectrum is still observed to change with 
wavelength. 

The KL region is quite massive and contains high 
column densities of both warm (T > 40 K) and cold 
(T < 25 K) dust. As a result, emission from cold and 
warm dust contribute nearly equally to the total flux ob- 
served at 350 and 450 /im. The coldest (T w 20 K) and 
warmest (T rs 45 K) regions of OMC-1 are regions to- 
wards KL and KHW (|Vaillancourtl l2002). As one moves 
towards KL from within 30" of the peak the tempera- 
ture of the cold component decreases by ~ 2 K, and 
the warm component increases by ~ 3K. These modest 
temperature changes are enough to shift the polarization 
minimum by as much as 70 - 130 /im (depending on the 
relative spectral indices of the components); moving to 
longer wavelengths towards the intensity peak. 

Therefore, we expect the polarization minimum to shift 
to longer wavelengths at the larger flux values towards 
the KL peak. If the observed wavelengths are near the 
minimum then we expect the 450/350 ratio to decrease 
from a value greater than unity when the minimum is be- 
low 400 /im, to a value approximately equal to one when 
the minimum is bracketed by the two observed wave- 
lengths, and to continue to decrease below one as the 
minimum moves to wavelengths greater than 400 fim. 
This is the trend observed in Figure 2] for points with 
P350 > 50%. 

4.3. Summary and Future Work 
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By comparing the 450/350 polarization ratio in 
OMC-1 with that in other bright Galactic clouds (Fig. [3]) 
we estimate that the minimum in the polarization spec- 
trum is not much less than 350 /im. However, given the 
variation in 60 - 1300 /im polarization spectrum from 
cloud-to-cloud, we do not draw strong conclusions about 
the polarization minimum in any cloud other than OMC- 
1. We are continuing our campaign to observe and an- 
alyze other bright clouds in a manner similar to that 
described here. 

The existence of a minimum in the far- 
infrared/submillimeter polarization spectrum, as 
opposed to a simple rise or fall from one end of the 
spectrum to the other, implies either (a) the existence 
of more than 2 dust temperature components, and/or 
(b) a change in the dust emissivity index as well as 
its temperature. Multi-wavelen gth observations in 
the far-infra r ed (5 - 200 ^m; IPotson et all 120001 : 
iDowell et ah! 120031 iVaillancourt et alj I2007D or sub- 



millimeter (> 450 am: ICurran &: Chrvsostomoul 120071 ; 
iBastien. Jenness. fe Molnarl 120051 : Matthews et al.. in 
preparation) at wavelengths on only one side of this 
minimum are insufficient for studying this behavior. 
SHARP (at 350 and 450 /im) therefore occupies a 
unique niche for studies of the polarization spectrum by 
providing a link between these two wavelength extremes. 
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AST 02-41356, AST 05-05230, and AST 05-05124. The 
CSO is supported by the NSF through grant AST 05- 
40882. 
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Fig. 1. — Polarimetric and photometric maps of OMC-1. Effective beamsizes (FWHM) for the photometric (9") and polarimetric (13") 
observations are shown as gray circles in (a) and (b). Coordinate offsets are measured with respect to OrilRC2 at 5 h 35 m 14.5 s , -5°22'31" 
(J2000). KL is the northernmost flux peak coincident with the coordinate origin and KHW/Orion-south is the peak ~ 1^5 to the south. 
Only polarization data satisfying P > 3a p are included, a) 350 (red) and 450 (blue) fim polarization vectors superposed on 350 /an flux 
contours. Contours are drawn at 2, 4, 6, 8, 10, 20, and 90% of the peak (« 780 Jy per 9" beam). 6) Inferred magnetic field vectors at 
350 and 450 /an drawn with a constant length (i.e., not proportional to the polarization amplitude); contours as in (a), c) The color scale 
shows the polarization ratio between the two wavelengths, P(450)/P(350). Contours at 350 fira are drawn at 4, 6, 10, 20, 30, 50, and 80% 
of the peak flux; the 50% contour is drawn thicker (see Fig.[4jl. 



